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ABSTRACT 



Rotation measurements in jets from T Tauri stars is a rather difficult task. 
Some jets seem to be rotating in a direction opposite to tliat of tlie underlying 
disk, although it is not yet clear if this affects the totahty or part of the outflows. 
On the other hand, Ulysses data also suggest that the solar wind may rotate in 
two opposite ways between the northern and southern hemispheres. We show 
that this result is not surprising as it may seem and that it emerges naturally 
from the ideal MHD equations. Specifically, counter rotating jets do not contra- 
dict the magneto-centrifugal driving of the flow nor prevent extraction of angular 
momentum from the disk. The demonstration of this result is shown by com- 
bining the ideal MHD equations for steady axisymmetric flows. Provided that 
the jet is decelerated below some given threshold beyond the Alfven surface, the 
flow will change its direction of rotation locally or globally. Counter-rotation is 
also possible for only some layers of the outflow at specific altitudes along the 
jet axis. We conclude that the counter rotation of winds or jets with respect 
to the source, star or disk, is not in contradiction with the magneto-centrifugal 
driving paradigm. This phenomenon may affect part of the outflow, either in 
one hemisphere, or only in some of the outflow layers. Prom a time dependent 
simulation, we illustrate this effect and show that it may not be permanent. 

Subject headings: Magnetohydrodynamics (MHD); Stars: pre-main sequence; Stars: 
winds, outflows; Stars: mass-loss; Stars: rotation; ISM: jets and outflows 
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Introduction 



Many a t temps have been made to measure jet rotation in various wavelengths (e.g. 



Coffey et al. 



20041 ). which is a rather delicate and difficult tas k. The RW Aur jet 



among other s , has been the most extensively studied such case. 



Cabrit et al. 



Coffey et al. 



(120061) and 



(120041 ) conclude that the rotation of the receding optical jet is opposite to that 



of the underlying disk. However, observations in the near UV do not confirm this result 

I II 1 

for the approaching jet (iCoffey et al.ll2012l) . This clearly shows the difficulties of rotation 
measurements which require a very precise experimental procedure. 

Even though the RW Aur jet is not a convincing case, counter rotation has been 
observed in several other jets. For inst ance, there is at lea st one knot (SKI) in HH212 that 



seems to counter rotate as explained in 



Coffey et al. 



( I2OIII ). with new evidence coming from 



Fe and HII lines. This is also the case for HHlll whose observations indicate a counter 
rotating knot too. 

Moreover, observations of the solar wind by Ulysses also show that in situ rotation 
measurements are delicate. 



Sauty et al. 



( 120051 ) have plotted in Fig. 1 of their paper, the 



latitudinal variation of the rotation velocity of the flow, the absolute value of which is very 
small. Although the conclusion is not definitive and could be due to some experimental 
uncertainties, the plot may also indicate that the solar wind in the northern hemisphere 
rotates in an opposite direction as compared to the wind in the southern hemisphere. 
Namely, the southern outflow may well be counter rotating with respect to the Sun. 

Despite the fact that rotation measurements in jets is a difficult task, counter rotation 
remains still an important open issue. First of all because according to the common 
view this may cont r adict the classical magneto-centrifugal outflow driving proposed in 



Blandford fc Paynd (119821 ) for jet launching from a Keplerian disk. The present letter 



aims at showing that this is not the case. Counter rotation is a direct consequence of the 
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velocity variation along the flow as well as of the flux tube geometry. It is shown that the 
conservation laws are satisfied and specifically that angular momentum flux is constant 
along the flow. This theoretical result is illustrated by a simulation which shows that non 
steady effects can also cause counter rotation. 



2. Governing equations for ideal steady axisymmetric MHD outflows 



2.1. Summary of the basic assumptions 



The basic equations governing ideal MHD outflows are the momentum, mass and 
magnetic flux conservation, the frozen-in law for infinite conducti vity and the firs t 



law of thermodynam ics. 



Heyvaerts &: Norman 



1982 



n steady and axisymmetric conditions (iTsinganosI 
19891 ) there are five integrals, quantities that are conserved along a 
given (poloidal) stream/field line. These are the magnetic flux A, the mass to magnetic flux 
ratio, the energy flux E{A), the angular momentum flux L{A) and the angular frequency 
or co-rotation frequency Q{A). In the following, we denote cylindrical coordinates with 
{z, zu, tp) and spherical coordinates with (r, 9, tp). 

The total angular momentum flux and the co-rotation frequency can be expressed in 
terms of the azimuthal and poloidal components of the magnetic and velocity fields. The 
toroidal component of the induction equation combined with the frozen flux condition gives 
the co-rotation frequency Q, 



Q{A) = -(v^- 

UJ \ Bp 



(1) 



The angular momentum flux is obtained by integrating the momentum equation in the 
azimuthal direction. 



L{A) = w{V, 



Br. 



ATxpVp 



(2) 
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The poloidal Alfven speed and the cyhndrical Alfven radius are defined as, 

For a smooth crossing of the Alfven surface, where Vp = Va, the cyhndrical radius along a 
given fieldline must adjust to this value, = ro*. 

From Eqs. ([I])-® we derive the azimuthal components as functions of the free integrals 
and the poloidal components of the fields, using the Alfven speed and radius: 

L W ro? L Vr, 'wi 



V = ^ ^ * y _ ^* IA\ 



With these expressions we show that a reversal of the rotation velocity is possible. Note 
that the toroidal components depend only on the poloidal ones and the free MHD integrals, 
being independent of the thermodynamics of the fiow. 



2.2. Reversal of the toroidal velocity 

Consider the fiow along a given stream/field-line as seen in Fig. [1] Assuming that the 

fiow remains everywhere super- Alfvenic after crossing the corresponding critical surface, the 

denominator of Eq. (jlj) is always positive, rryj — 1 > 0. The sign of the toroidal velocity is 

Va 

then given by the numerator, 

'V^ w^' 

sgn(V;) = sgn ( 2 1 • (5) 

Va ^* 

The initially positive rotation velocity will change sign if the velocity of the fiow along 
the fiux tube is less than the threshold value, 

< Va = threshold ■ (6) 
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Fig. 1. — 3D view of a flux tube along the flow, where the cylindrical radius is shown at the 
Alfven transition and another point downstream. On the right part of the figure, a zoom of 
the poloidal plane shows the definition of V^, 5S, and ^-S'^;. 
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Note that this value is always above the local Alfven speed because we expect > ro^. 
Along a given flux tube, we can write the magnetic and mass flux conservation as 

BpSS = B^6S-i, = constant (7) 

and 

pVpSS = Pi,Vi,6Si, = constant . (8) 

Note that 6S is the surface element perpendicular to the poloidal velocity (see Fig. [1]), such 
that one can always replace Vp by any of its components using the following equation, 

Vp-6S = Vp6S = VJS, = VrSSr = V,6Si , (9) 

where Vi is the i-th component of the velocity and dSi is the cross section perpendicular 
to the z-th direction. The calculation is simpler if one works with the poloidal velocity. 
Nevertheless, it is more straightforward to compare the z component of the velocity and 
the cross section 63^, which is in the zu direction, since these are the observed quantities. 

Using these two equalities in the numerator of the azimuthal velocity we successively 

get: 



T/2 2 ~ r>2 2 ~ rc2 „2 

^* ^* ^* 

2 _ d2 



(10) 



Since at the Alfven surface we have Anp^^Vj = B^, the above expression simplifies further 

T/2 2 ^^* V 2 XC T/ ^2 



Thus, the sign of this last expression determines the sign of the azimuthal velocity in the 
super Alfvenic region: 



sgn(l/.)-gn(^-^^)^ (12) 
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This gives a new expression for the threshold value defined in Eq. ([6]), 

^threshold = K — 2T7T ■ (13) 

If the fiow remains super Alfvenic, reversal of the rotation takes place when the velocity 
drops below some threshold value. The fiow can decelerate either because it expands and 
this provokes adiabatic cooling, or because its kinetic energy drops due to some other 
mechanism such as radiative losses. The reversal might also occur if the threshold value 
increases. This may be the case when the cross section of the tube decreases significantly 
while the fiux tube itself widens. This is illustrated in Fig. [2l namely, the toroidal velocity 
becomes negative at the location where the poloidal speed curve crosses the line of the 
threshold velocity. 

Counter rotation is also obviously obtained if the fiow becomes sub-Alfvenic. In this 
case, the denominator of Eq. (jl]) becomes negative and the reversal of the toroidal velocity 



is accompanied by a reversal of the toroidal magneti c field too. 



Fendt 



'his may happen after a 



( 120111 ). that the reversal of the 



shock, which is probably similar to the suggestion of 
toroidal field is due to the presence of shocks. 

Last, we also note that in the sub-Alfvenic part of the fiow close to the source, rotation 
may also change sign if there is a strong local acceleration and if the velocity increases 
above the threshold value (instead of below in the super- Alfvenic part). 

Another way of exploring the counter rotation possibilities is to express the numerator 
of the rotation speed from Eq. (jlj) using the density ratio, 

(14) 

2 

The sign of this quantity determines the sign of V"^, i.e. reversal will occur if p > p*^. 



Fig. 2. — Sketch along a given streamline showing the evolution of various velocities: the 
poloidal component V^, the threshold value T^hreshoid; the poloidal Alfven speed Va and the 
rotation speed V^. The threshold value varies due to geometrical effects but remains close to 
\4 (dotted line), which is the value of the poloidal speed at the Aflven point. The rotation 
changes direction at the location where the poloidal velocity far from the launching region 
becomes lower than the threshold value. Note that the flow remains super- Alfvenic after the 
Alfven point. 
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2.3. Discussion on the threshold value 

For meridionally self-similar models, it is more convenient to use the radial component 
of the velocity, i.e. Eq. iQ. The criterion for a change of sign in the flow rotational velocity 
is now given by, 

cos 6 

Vr < K,. ^ . 15 

This is particularly interesting for analytical stellar jet or wind models. 

Similarly for radially self-similar models, reversal of the rotation occurs where, 

|^.|<|V^.,.|^. (16) 
sm 

This formula applies instead on analytical disk wind models. 

If the physical quantities are uniform across the jet then SS is proportional to ro^. The 
condition then reduces to the much simpler expression, V < V^. The threshold value is 
exactly the Alfven speed at the Alfven surface. 

To compare this criterion with observations and numerical simulations it is useful 
to work with the vertical component of the velocity using Eq. (|9]) because this is the jet 
velocity component measured. The surface element is then replaced by its projection 
along z, dSz, which is the cross section of the flux tube in the horizontal tz7 direction that 
corresponds to the jet radius as it is measured, 

v;<v;4§^. (17) 

If we average the quantities at a given altitude z, we may have a criterion for the 
averaged velocity using the previous equation, namely. 



v. = '- — —<^ — — = v.. [li 
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The criterion (i.e. V< V^) is equivalent to the criterion for transversaly uniform flow. 



This may be more adequate to analyse observational dat a and jet models such 



as the one proposed by 



Lovelace. Berk, fc ContopoulosI (Il99ll ). In fact, Eq. (16) of 



Lovelace. Berk, fc ContopoulosI ( 1l99ll ) shows that an increase of the average jet velocity 
induces a decrease of the rotational one, which could become negative. 



2.4. Counter rotation in simulations 



We have performed a numerical simulation in which several reversals of the toroidal 



Matsakos et al. 



torn 



velocity are observed. We have adopted the initial conditions of 
model NIPH, in which two analytical solutions for the stellar and disk wind components 
are combined and a constant external pressure is applied beyond a speciflc radius. We have 



carried out the temporal integration of the MHD equations using PLUTO ( iMignone et al. 



20071 ). in a resolution of 512 x 2048 up to a flnal time that corresponds to ~100 crossing 



times of the simulation box. 

During the initial transient phase, a strong shock propagates upwards due to the fact 
that the initial state is not in equilibrium. Together with the shock, a signiflcant reversal 
of the toroidal velocity is also observed, although this phenomenon is not directly related 
with the above analysis. Once the system reaches a quasi-steady state, several layers 
demonstrate an opposite V^. These counter rotating regions remain roughly at the same 
location for large time intervals, before they vanish and reappear. They obviously remain 
transient phenomena. 

The left panel of Fig. [3] displays contours of the toroidal velocity on the left and the 
poloidal component on the right. We have selected a streamline {-cu{z), grey and white 
solid lines of left panel) that crosses the region where the rotation reverses, along which the 
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Fig. 3. — Snapshot of a two component jet simulation surrounded by an external medium 
of constant pressure. On the left panel, the toroidal velocity (left) and the corresponding 
poloidal velocity (right) are displayed. It also gives the cylindrical radius ix>{z) as a function 
of z. The reversal of the rotation corresponds to the darkest zone of the contours. A 
specific streamline has been selected (solid line) along which the quantities shown on the right 
panel are plotted. From top to bottom: the cylindrical radius ixj{s) the angular momentum 
flux, the toroidal speed, the poloidal and Alfven speeds, and the poloidal Alfvenic Mach 
number as functions of the curvilinear abscissa, s. The angular momentum flux (solid line) 
has two components, the hydrodynamical one (I^hydroj dotted line) and the magnetic one 
(-^MAG) dashed line). 

cylindrical radius, the angular momentum fluxes, the toroidal, poloidal and poloidal Alfven 
speeds as well as the poloidal Alfvenic Mach number are shown on the right panel of Fig. El 
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The plots suggest that the component becomes negative at the location where 
the poloidal velocity decreases approaching the Alfven speed. Note that the flow remains 
always super-Alfvenic. The dotted dashed vertical lines indicate the change of sign of the 
rotation. In this region, the poloidal velocity reaches the local threshold value which varies 
along the flow. Although the jet is not an exact steady flow, the simulation results can be 
explained based on the analysis that describes the super Alfvenic regime. 

On Fig. [3] we plot the total angular momentum flux L and its two components 
(see Eq. [2]), namely th hydrodynamical part Lhydro = '^V^p and the magnetic one 
-^MAG = ~^ 4^V '^^^ total angular momentum flux is clearly not conserved because 
the configuration is not in an exact steady state. However the reversal of the toroidal 
velocity corresponds to an exchange between the two components where the magnetic part 



becomes dominant. L drops rapidly up to s = 190. In this first part the reversal o: 



Fendt 



the 



fl201lh 



velocity is strongly time dependent and related to the shock as suggested by 
Above s = 190, L is almost constant and in this region our steady state criterion for the 
reversal is valid. It is a nice illustration of the above theory. 



3. Conclusions 

We have presented a simple study showing that counter rotation may occur in any 
jet or wind independently of the launching mechanism. Counter rotation may be observed 
in the outer layers of a disk wind, as shown in our numerical simulation and is observed 
in some protostellar jets. It may also occur in a stellar wind, as suggested from the data 
of the solar wind. This effect does not contradict the disk launching mechanism, nor the 
conservation of the angular momentum along the poloidal field lines. The criterion that 
marks the reversal of is determined by either a deceleration of the fiow or a specific 
geometrical configuration. 
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For the Sun, the toroidal velocity of the solar wind is very small or negative because 
there is almost no acceleration after the Alfven point. In outflows associated with YSOs we 
expect that there may be reversals of the rotation in some of the layers of the flow but not 
necessarily across or along the whole jet. There is no evidence that this phenomenon should 
be permanent even if it is not a transient modification of the flow dynamics. In all cases, 
the total angular momentum is conserved, Eq. ([2]), although the kinetic and magnetic terms 
constituting L{A) are expected to exchange the amount they store. The situation is similar 
to the transfer of energy from the enthalpy to the kinetic term in the Bernoulli integral in 
the thermally driven Parker wind, or, to the transfer of energy from the Poynting energy 
term at the base to the kinetic energy term at large distances in the generalized Bernoulli 
integral in magneto-centrifugal winds. 
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